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RADIATION  FROM  INTENSE  ELECTRON  BEAMS  ASSOCIATED  WITH 
THE  CERENKOV  MECHANISM 

1.  INTRODUCTION 

Normally  a  charge  or  system  of  charges  moving  at  constant 
velocity  does  not  radiate  electromagnetic  energy,  but  if  the 
velocity  of  the  charge  is  greater  than  the  velocity  of  light  in  a 
medium,  radiation  occurs.  It  is  broad  band;  the  intensity 
increases  with  frequency;  and  the  propagation  is  confined  to  a 
cone,  much  like  a  shock  wave  in  air.  Equations  (?)  and (3)  in  a 
later  section  describe  these  results.  In  previous  work  l"4,  we 
described  theory  and  experiments  showing  that  this  mechanism 
produces  microwaves  from  bunched  electron  beams,  in  addition  to 
the  usual  optical  radiation,  as  studied  by  Cerenkov  and  many 
others . 

In  this  report  we  consider  aspects  of  coherent  Cerenkov 
radiation,  with  emphasis  on  radio  and  microwave  frequencies.  As 
will  be  shown,  power  levels  may  become  very  high,  in  the  megawatt 
range,  so  that  serious  consideration  should  be  given  to  this 
mechanism.  The  scopjp  of  this  report  is  to  outline  the  systematics 
of  a  theory  in  which  the  dielectric  (which  lowers  the  velocity  of 
light)  is  considered  to  be  linear  and  non-dispersive .  The  aspects 
considered  which  are  not  treated  in  the  extensive  literature  are 
calculations  for  electrons  occurring  in  finite  size  bunches,  both 
single  and  repeated.  Future  work  should  aim  at  the  effect  of 
ionization  in  the  dielectric  medium,  betatron  oscillations  of  the 
beam  electrons,  and  other  mechanisms  which  could  seriously  modify 
or  degrade  the  radiation. 


2.  SUMMARY  OF  PREVIOUS  V7QRK 

In  the  earlier  work1-4,  it  was  assumed  that  electrons  were 
emitted  from  an  accelerator  in  bunches,  described  by  a  charge 
density  p{r),  and  all  had  the  same  velocity  v  such  that  the  bunch 
did  not  change  shape  in  time.  The  current  density  is  then  given 
by 

J(r)  =  p  v  (1) 


Eq.  (2)  refers  to  a  single  bunch,  the  radiation  has  a 
continuous  frequency  spectrum,  and  the  formula  gives  the  energy 
radiated  at  angular  frequency  oj  in  a  range  dco,  per  unit  length. 
Eq.  (3)  is  the  energy-per-unit  path  per  bunch  for  the  case  of 
periodic  bunches,  emitted  at  an  angular  frequency  u)0.  Here  u>  is 
a  harmonic  of  uiq,  in  contrast  to  the  continuous  frequency 
spectrum  of  Eq.  (2).  In  both  cases,  q  is  the  total  charge  of  one 
bunch,  and  F(k)  is  a  form  factor  for  the  charge  distribution 

FOt )  =  q'1///  p(r)  e^*r  dV  (4) 


9C  is  the  Cerenkov  angle, 


q  is  the  charge  per  bunch,  and  k  is  propagation  vector  for  the 
emitted  radiation,  such  that  k  =  u/c,  and  the  direction  is  along 
the  Cerenkov  cone . 

B.  Polarization 

The  electron  field  is  in  the  plane  of  the  observer  and  the 
electron  path,  as  it  is  for  usual  optical  Cerenkov  radiation. 

C.  Coherence 

If  the  charge  is  a  point  charge,  F  -  1  and  in  fact,  for  a 
finite  bunch,  we  will  have  F  ~  1  for  k£  <<  1,  where  k  *>  2n/X. 

This  corresponds  to  X  >  l,  where  i  is  the  size  parameter  of  the 
bunch.  For  these  long  wave  lengths,  either  formula,  Eq.  (2)  or 
(3),  shows  that  the  intensity  is  found  by  setting  q  =  ne  where  n 
is  the  number  of  electrons  in  the  bunch.  Thus  the  intensity  is 
proportional  to  n2e2.  As  higher  frequencies  are  examined,  the 
factor  F  becomes  small,  and  finally,  for  very  high  frequencies 
(optical  range),  we  would  expect  incoherent  radiation,  in  which  q2 
is  replaced  by  ne2. 

D.  Diffraction  Effect 

An  important  result  is  that  if  the  region  of  emission  of 
radiation  is  finite,  of  length  L,  the  emitted  radiation  is  no 
longer  confined  to  the  Cerenkov  cone  9  *  0C  but  is  spread  about 
that  angle  and  the  spreading  depends  on  the  length  of  the  path  in 
the  gas  and  the  wave  length  of  the  radiation.  It  is  more  closely 
related  to  the  diffraction  pattern  of  an  end  fire  antenna,  rather 
than  the  broadside  array;  the  latter  is  related  to  the  single 
slit  of  optics.  This  effect  is  independent  of  bunch  size  and 
repetition  rate.  Associated  with  diffraction  is  an  enhancement, 
qualitatively  described  below. 
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E. 


Enhancement  from  Diffraction 


In  most  diffraction  or  interference  situations  in  optics, 
the  energy  is  usually  redistributed  in  angle  if  parameters  are 
changed.  For  example,  if  the  width  of  a  single  slit  is  changed 
but  the  total  power  passing  through  the  slit  remains  constant,  the 
total  power  integrated  over  the  diffraction  pattern  on  a  screen 
remains  constant.  In  the  Cerenkov  case,  the  development  of  Ref.  3 
showed  that  for  small  Cerenkov  anqles,  the  intensity  at  0C 
remains  constant  as  L  decreases,  but  the  total  power  radiated  at 
angles  0  >  9C  increases  dramatically.  The  total  power  radiated 
increases  by  factors  of  up  to  50  in  our  experimental  situation. 

F.  Form  Factor  Effect 

The  basic  Cerenkov  radiation  strength  tends  to  increase 
with  frequency,  for  either  coherent  or  incoherent  emission. 
Starting  at  low  freauencies,  the  coherent,  emission  increases  with 
w  until  F  becomes  small,  at  u  *  u^.  At  this  frequency,  the  wave 
length  is  short  enough  so  that  radiation  emitted  from  electrons  at 
the  front  and  rear  of  the  bunch  begin  to  interfere  destructively. 
This  leads  to  destructive  interference  and  lowered  emissions  of 
power  above  <%.  Experiments  in  principle  could  demonstrate  this 
effect,  but  at  this  time  we  have  been  unable  to  make  measurements 
in  the  high  frequency  microwave  ranoe  ( \  ~  2  mm)  which  would 
quantitatively  demonstrate  the  predictions. 

As  a  further  remark,  it  should  be  possible  with  precise 
measurements  to  determine  size  parameters  of  the  electron  bunches 
at  a  distance,  or,  in  other  words,  to  do  beam  diagnostics  by  this 


method.  These  possibilities  will  be  explored  in  another  report. 
3.  Approximate  Results  for  Radiated  Power 


i 


The  results  mentioned  previously  may  be  related  to  more 
conventional  parameters  of  pulse  length  and  peak  current.  If  we 
assume  a  pulse  in  the  form  of  a  uniform  cylinder,  J  =  pv  and  the 
current  I  =  JA,  q  =  pAi  =  pAvT  where  l  =  the  length  of  the  bunch 
and  T  *  the  bunch  length  measured  as  a  time.  Thus  we  find 

q  -  IT  (6) 

For  a  Gaussian  pulse,  the  charge  density  is 

-(x2  +  y2)/a2  -b  2z  2/b  2  (7) 

p  —  p  q  e  e 

It  may  be  shown  that  the  peak  current  in  the  middle  of  the 
pulse  (at  z  =  0)  is 


I 

I 


I  = 
o 


irb 


and  that  the  normalization  constant 


charge  of  the  bunch  by 


q  =  p_ 


3/2  2 

it  a  b 


(8) 

PQ  is  related  to  the  total 

(9) 


Then  we  may  let  the  effective  length  of  the  bunch  be  Z  =  bn 

and  use  q  =  I0T  =  I0  Jl/v,  just  as  for  the  uniform  cylinder. 

Now  consider  the  sum  over  frequencies  in  Eq.  (3).  Assuming 

that  F(k)  is  constant  up  to  some  k  which  corresponds  to 

c  max 

i»)  =  Nu)n,  we  find 

m  0 
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J>U)0  =  Ino)02 

=  N(N+1)  u>02/2 
*  N2uj02/2 
=  <%2/2 

Thus 


E““o  =  uhi2/2 


(10) 


For  a  single  bunch,  the  sum  over  the  continuous  frequency 
distribution  (2)  gives 

/  mu)dto  =  a)  2/2  (11) 

6  m 

Thus  the  results  of  Eqs.  (10)  and  (11)  above  mav  be  used  for 
either  single  or  repeated  pulses  to  calculate  the  total  power 
radiated  at  all  frequencies. 

Now  we  may  use  the  results  of  Eq.  (6)-(t?)  to  re-express  the 
results  of  Eq.  (?)  or  (3).  The  latter  may  be  multiplied  by  v  to 


convert  to  power,  then  summed  or  integrated  over  w  for  the  total 


power  emitted,  which  is 

trr  „  ,  2  2  2 

P  =  -JLo.  v  s i n 2 0  In  T  f 
2  c  0  m 

This  holds  for  either  periodic  or  a  single  bunch. 


(12) 

In  either 


case  it  is  the  rate  at  which  electromagnetic  radiant  energy  is 


emitted  while  the  bunch  is  in  flight,  per  bunch  (in  the  case  of 
periodic  bunches). 

Noting  that  fm  =  1/T  holds,  which  is  equivalent  to  the 
uncertainty  relation  in  quantum  theory  or  the  bandwidth-pulse 
length  relation  in  communications. 
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•  1 


K.J 


P  =  JLHj 


v  sin20c  I02 


Thus  the  power  emitted  for  coherent  radiation  is  proportional 
to  the  square  of  the  beam  current.  This  makes  explicit  one  of  the 
main  points  of  this  report:  that  the  emitted  power  increases 
rapidly  with  the  beam  current.  Thus,  although  weak  beams  emit 
little  radiation,  the  radiation  becomes  important  for  large 
currents.  When  we  note  that  the  radiation  emitted  is  proportional 
to  I02  but  the  energy  in  the  beam  is  proportional  to  IQ,  so  that 
relative  efficiency  as  a  radiator  increases  in  proportion  to  IQ. 
This  may  be  made  quantitative  by  multiplying  (13)  by  L/v  to  oet 
the  total  energy  emitted  by  the  bunch  during  its  flight  through  a 
path  L,  where  L  is  the  range.  The  energy  put  into  the  bunch  is 
I g T ,  so  that  the  efficiency  e,  defined  as  the  ratio  of  the  energy 
radiated  (Cerenkov)  to  the  total  energy  of  the  beam  becomes 

e  =  £$I_  sin\  -I*  (14) 

where  VQ  is  the  energy  of  the  beam  expressed  in  volts  and  L  is  the 
range  of  the  beam  in  air.  Strictly  speakinq,  L  is  the  distance  at 
which  Cerenkov  radiation  ceases.  The  precise  value  of  L  would 
depend  on  the  relative  importance  of  bremsstrahlung,  ionization, 
and  instabilities  which  could  disperse  the  electron  bunch,  and 
these  aspects  should  be  considered  in  greater  detail . 

The  overall  systematics  of  the  radiation  associated  with  the 
Cerenkov  mechanism  are  expressed  by  Eq.  (12),  (13),  and  (14).  The 

following  features  should  be  noted: 
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A.  The  power  radiated  per  bunch  (eq.  12,  13)  is  proportional 
to  the  square  of  the  peak  beam  current,  I 02,  and  is  independent  of 
the  pulse  length  T.  This  depends  on  the  approximation  that 

fmT  =  1,  where  T  is  the  pulse  length  and  fm  is  the  maximum 
frequency  for  which  the  fourier  transform  of  the  pulse  has 
significant  strength.  This  expression  for  the  power  radiated  is 
the  energy -per-unit  time  while  the  bunch  is  in  flight,  from  the 
accelerator  to  the  point  where  it  ceases  to  radiate  effectively. 
The  latter  point  in  the  simplest  consideration  is  the  point  at 
which  the  bunch  has  lost  enough  energy  so  that  it  is  below  the 
threshold  for  Cerenkov  radiation,  which  is  20  MeV  for  electrons  in 
normal  air. 

B.  As  T  is  decreased,  assuming  IQ  is  constant,  the  power 
radiated  remains  constant.  But  the  frequency  at  which  most  of  the 
radiation  occurs  is  increased  as  T  is  decreased. 

C .  The  efficiency  ( Eq .  14)  depends  on  Ic/T  showing  that  for 
high  currents  and  short  bunches,  the  Cerenkov  mechanism  could  not 
only  lead  to  sizable  amounts  of  radiation  but  also  dominate  the 
other  modes  of  energy  loss. 

D.  As  a  guide  to  what  may  be  expected  in  the  case  of  plasma 
shielding,  we  might  expect  that  some  length  l'  of  the  nose  of  the 
beam  bunch  might  be  unshielded  by  the  plasma.  From  the  above 
considerations,  the  lenqth  l'  would  not  change  the  total  power 
radiated,  but  would  change  dramatically  the  freauency  at  which  the 
radiation  occurs.  This  is  because  if  only  the  nose  region  of 
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of  length  V  radiated,  the  total  effective  radiating  charge  would 
be  proportional  to  But  the  frequency  of  the  maximum  radiation 

would  increase  as  1 / 1‘ ,  so  the  power  radiated  would  remain 
constant  as  Z‘  is  varied. 
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APPENDIX  A  -  EXTREME  VALUES  OF  THE  RADIATION  FUNCTION 


The  radiated  power  per  unit  solid  angle  W(v,n)  can  be  written  as  the 
product  of  two  factors: 

W  =  Q  R2  (Al) 

where  Q  is  a  function  having  the  dimensions  of  power  (watts)  and  R  is  a 
dimensionless  radiation  function. 

**  V  2 

Q - q2  ( A2) 

8tt2 

R  *  (kL)  SinO  F(k)  I(u)  (A3) 

Here,  p  is  the  permeability  of  the  medium,  c  is  the  velocity  of  light  in 
the  medium,  and  q  is  the  charge  of  an  individual  electron  bunch.  If  an  RF 
Linac  is  the  source  of  the  electrons,  the  frequency  v0  is  that  of  the  Linac ; 
however,  in  all  cases,  vQ  corresponds  to  Che  spatial  periodicity,  l  *  c/v0, 
of  the  electron  bunches. 

In  (3),  k  is  the  wave  vector  of  the  emitted  radiation,  L  is  the  beam 
interaction  length,  and  F(k)  is  the  form  factor  defined  such  that 


g0’(k)  =■  q  F(k) 


(A4) 


and  £g'(k)  is  the  Fourier  transform  of  the  charge  distribution  of  an 

charge  bunch.  The  diffraction  function 
j.  m  Sin  u 
u 

is  a  function  of  the  diffraction  variable 

u  =  kL  (Cos9  -  CosOl 

~  c 


individual 

(A5) 

(A6) 


and  consequently,  R  is  a  somewhat  comolicated  function  of  0  and  kL. 
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MINIMUM  VALUES 


By  inspection,  it  is  obvious  that  R2  and  thus  W  will  have  minima  when 
I(u)  *  0.  These  true  zeroes  of  the  function  occur  when  u  is  equal  to  integral 
values  of  it;  i.e.,  u  =  mu.  Substituting  this  condition  in  (A6)  leads  to  the 

angular  values  for  the  zeroes  in  the  radiated  power 

Cos9  =  Cos 3  -  m  '  (A7) 

c  VcL 

which  can  be  written  in  terms  of  the  wavelength  (in  the  medium)  of  the 
radiation  as 

CosO  =  Cos 0  -  m  (^)  (A8) 

This  relation  is  exact  and  gives  the  zeroes  of  R  for  all  allowed  values  of 
m  for  which  Cos0  <  1,  I(u)  is  a  maximum  at  m  =  0  (which  gives  u  =  0),  and  the 

maximum  in  R  tends  to  occur  somewhere  in  that  vicinity.  The  first  zero  is  that 

one  which  is  adjacent  to  the  maximum  and  is  given  by  (A8)  when  m  *  +1.  Values 
of  this  angle  as  a  function  of  y  are  shown  as  a  Log-Log  plot  in  Fig.  1.  The 
relationship  between  0  and  is  approximately  that  of  the  inverse  square  root 
over  the  range  10  <  L/  A  <  100.  Within  this  range,  the  angular  value  (in 
degrees)  of  the  first  minimum  is  given  by  the  empirical  relation 

9l  =  -,8ri.  ( A9) 

VL/A 

The  physical  range  of  0  is  from  0  to  tt.  Substitution  of  the  successively 
larger  positive  values  of  m  =  2,  3,  4...  into  (A7)  or  (A8)  gives  a  set  of 
increasing  angles  for  the  functional  zeroes.  If  A/L  «  1,  many  zeroes  will 
occur  before  the  physical  limit  is  reached. 

In  addition,  if  A/L  «  1,  there  is  the  possibility  of  having  functional 
zeroes  at  angular  values  less  than  that  of  the  main  peak.  These  occur  with  the 
substitution  of  successive  negative  m  values  into  (A7)  or  (A8) ,  giving  rise  to 
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subsidiary  radiation  maxima  between  the  beam  line  and  the  main  diffraction 
peak. 

MAXIMUM  VALUES 

To  find  the  maximum  values  of  ^  one  proceeds  in  the  usual  fashion  by 
taking  the  derivative  and  setting  it  to  zero. 


3W 


3R 


“  2  Q  R  To  '  0 

This  equation  is  satisfied  for  either  R  *  0  or 


3R 

33 


(A10) 

0.  The  first  condition  is 


that  considered  above  for  the  zeroes  of  the  function.  The  remaining  one, 

3R 

*  0,  pertains  to  the  maxima.  Thus,  it  is  sufficient  to  consider  maximum 
values  of  R  which  will  also  be  maximum  values  for  W  ~  R2. 

From  (A3),  the  derivative  of  R  is 

-S  *  kL  {  sin 9  [f  4i  ♦  I  -«]  +  FI  Cos0  }  =  0  (All) 


and  the  derivative  of  the  diffraction  function  I(u)  can  be  written 


31  dl  3u  kL  c.  _ 

is ’at  is"  1  —  s,n0 


(A12) 


Substituting  this  value  in  (All)  and  discounting  the  solution  for  kL  =  0  gives 


the  condition  for  maxima 


{  Sin9  [I'  F  ~  Sin9  +  I  ~]  +  FI  Cos©}  =  0 


3F  ■ 


3Q- 


(A13) 


In  order  to  proceed  further,  the  functional  dependence  of  the  form  factor 


F  must  be  known. 


A  gausian  form  factor  has  been  assumed  in  this  and  other  reports. 

F(k)  “  Exp  (-a2  Cos  20)  (A14) 

where  a  *  kb/4  corresponding  to  a  one-dimensional  charge  distribution 
>  Z2 

Po'  (r)  =  q  Exp  ( - )  (A15) 


The  derivative  of  this  form  factor  is 

H  -  si» 23  p 


(A16) 


vrfiich  is  small  for  the  parameters  of  interest  in  Linac  experiments;  i.e., 
b  *  0.24  cm,  microwave  frequencies,  and  angles  up  to  approximately  25*. 
Therefore,  a  first  approximation  to  solving  (A13)  can  be  obtained  by  setting 
the  derivative  of  the  form  factor  in  (A13)  to  zero.  The  resulting  approximate 
maxima  condition  is  independent  of  F  and  can  be  rearranged  so  that  I  and  its 


derivative  are  on  the  same  side. 

-  r'<u>  ,  !_  Cos  3 
l(u)  kL  Sin 2  9 

The  ratio  of  I'  to  I  can  be  expanded  as  a  power  series  in  u  as 


(Al?) 


3  45 


(A18) 


which  converges  rapidly  for  small  values  of  u.  Substituting  the  first  term  in 
the  expansion  into  (A17)  and  subsequently  eliminating  the  diffraction  variable 
through  its  definition  (A6)  leads  to  an  approximate  equation  for  the  maxima, 


which  is  cubic  in  Y  *  CosO. 

-Y3  “  Cos 0  Y2  -  [l  ♦  —  ]  Y  +  CosO-  *  0 


(A19) 


In  this  equation,  the  beam  parameters  enter  through  the  charge  velocity  in 
Cos  0C  *  c/v  and  the  radiation  wavelengths  in  x  8  kL  *  2n  L/X.  Numerical 
solutions  can  be  carried  out  on  a  hand  calculator  with  results  within 


approximately  0.3’  of  those  obtained  from  large-scale  computer  calculations. 
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Results  for  the  major  radiation  lobe  (the  first  maximum)  are  also  plotted  in 
Fig.  1.  These  also  show  a  straight  line  dependence  on  a  Log-Log  plot  and 
therefore  can  also  be  represented  as  a  power  law 


which  is  accurate  to  approximately  5%. 


APPENDIX  B  -  EFFECT  OF  BEAM  LENGTH  ON  THE  RADIATION  FUNCTION 


In  previous  reports,  several  calculated  curves  of  a  radiation  function 
which  differs  only  by  constants  from  R2  were  presented.  Fig.  2  is 
reproduced  from  Fig.  5  of  Reference  4  and  relabeled  to  be  consistent  with  the 
notation  of  Appendix  A.  Fig.  2  shows  the  vartiation  of  the  radiation  function 
for  8.55  GHz  radiation  from  the  NPS  Linac  for  beam  interaction  lengths  in  the 
vicinity  of  100  cm.  As  noted  previously  in  this  and  other  reprots,  increasing 
L  tends  to  move  the  main  diffraction  lobe  toward  the  beam  line,  and  Appendix  A 
gives  an  approximate  formula  for  calculating  the  position  of  the  first  lobe. 
Fig.  3  (to  different  scales)  shows  this  effect  for  larger  beam  paths  from  an  S- 
Band  Linac.  As  L  increases,  the  principle  lobe  moves  in  toward  the  Cerenkov 
angle  9C;  i.e.,  ^  >  3C,  and  the  radiated  power  per  unit  solid  angle  at 
the  maximum  increases  greatly. 


17 


FIGURE  CAPTIONS 


Fig.  1  Cerenkov  radiation  extremum  angles  in  degrees  vs  L/A.  The  wavelength  A 
of  the  Cerenkov  radiation  is  measured  in  the  medium  in  which  the  path 
length  is  L.  Curve  B  is  for  the  first  minimum,  the  solution  to  (A8) 
with  m  *  +1.  Curve  A  is  for  the  first  maximum,  the  solution  to  the 
approximate  equation  (A19).  Over  the  range  shown,  A  is  straight  and  B 
is  slightly  curved.  Equations  for  the  curves  are  given  in  the  text. 

Fig.  2  Calculated  third  harmonic  (8.55  GHz)  radiation  patterns  produced  by 

electron  bunches  from  an  S-band  (2.85  GHz)  Linac.  The  path  lengths  in 
air  are  70,  90,  and  150  cm.  Vertical  scale  is  arbitrary  and  the  same 
for  all  three  curves. 

Fig.  3  Calculated  third  harmonic  radiation  patterns  produced  by  electron 

bunches  from  an  S-band  Linac.  The  solid  curve  is  for  a  path  in  air  of 
1000  cm,  and  the  dotted  curve  is  for  a  path  length  of  10,000  cm. 
Vertical  scale  is  arbitrary  and  the  same  for  both  curves  but  much 
reduced  from  that  of  Fig.  2. 
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